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©Modeling Electromagnetic Scattering From Complex Inhomo- 
geneous 0 bjects 

Complex, inhomogeneous objects can be easily modeled using commercial CAD packages. 

NASA Langley Research Center, Hampton, Virginia 


This software innovation is designed 
to develop a mathematical formulation 
to estimate the electromagnetic scatter- 
ing characteristics of complex, inhomo- 
geneous objects using the finite-ele- 
ment-method (FEM) and 
method-of-moments (MoM) concepts, 
as well as to develop a FORTRAN code 
called FEMOM3DS (Finite Element 
Method and Method of Moments for 3- 
Dimensional Scattering), which will im- 
plement the steps that are described in 
the mathematical formulation. Very 
complex objects can be easily modeled, 
and the operator of the code is not re- 
quired to know the details of electro- 
magnetic theory to study electromag- 
netic scattering. 

Consider a complex inhomogeneous 
object geometry, the electromagnetic 
scattering characteristics of which are 
to be estimated when illuminated by a 
plane electromagnetic wave. To facili- 
tate the mathematical formulation 
using the hybrid FEM / MoM , the object 
is assumed to be enclosed by a surface 
indicated as a fictitious outer bound- 
ary. Inside the boundary, the electro- 
magnetic fields are obtained using the 


FEM, whereas the electromagnetic 
fields outside the fictitious boundary 
are obtained using the assumed equiva- 
lent electric and magnetic currents 
flowing on the fictitious boundary. 
Continuity of the electromagnetic 
fields across the fictitious boundary re- 
sults in partly sparse/ dense matrix that 
is solved for the unknown fields inside 
the fictitious boundary including the 
boundary surface. 

The steps laid out in the mathemati- 
cal formulation are carried out in 
FEMOM3DS. Along with the main 
FEMOM3DS, the innovation uses a 
commercial computer-aided-design 
(CAD) package for geometrical model- 
ing, and a post-processing package such 
as Techplot for displaying graphically 
the results obtained using FEMOM3DS. 
In the present FEMOM3DS code the 
COSMOS/ M commercial software is 
used as a CAD tool to model the geom- 
etry of a given problem. The COS- 
MOS/M is also used to discretize the 
FEM region using the tetrahedral ele- 
ments. Various inhomogeneous regions 
are taken care of by having many parts 
in the FEM region. Using the COS- 


MOS/M, common boundaries where 
boundary conditions are to be imple- 
mented are also identified. The data 
file created by the COSMOS/ M for 
node and element information is then 
generated and processed through the 
preprocessor part of FEMOM3DS code 
to create edge information and node 
information. The preprocessed data are 
then run through the main part of 
FEMOM3DS to obtain electromagnetic 
scattering. The output files from the 
FEMOM3DS can be used for displaying 
the results in a graphical format. 

The FEMOM3DS is written in FOR- 
TRAN 77. The code was successfully 
compiled on a CONVEX machine. 
Flowever, the code can be compiled on 
any 32-bit machine like PCs or SUN, 
SGI UNIX Station. To get correct re- 
sults, dimensions must be given in cen- 
timeters, frequency of operation must 
be given in GFIz and incident and ob- 
servation angles must be specified in 
degrees. 

This work was done by Manohar Desh- 
pande and C. J. Reddy of Langley Research 
Center. Further information is contained in a 
TSP (see page 1).LAR -17 090-1 


©Visual 0 bject Recognition and Tracking of Tools 

This method can be used to track tools held and used by humans, such as surgical tools. 

Lyndon B. Johnson SpaceCenter, Houston, Texas 


A method has been created to auto- 
matically build an algorithm off-line, 
using computer-aided design (CAD) 
models, and to apply this at runtime. 
The object type is discriminated, and 
the position and orientation are identi- 
fied. This system can work with a single 
image and can provide improved per- 
formance using multiple images pro- 
vided from videos. 

The spatial processing unit uses three 
stages: (1) segmentation: (2) initial type, 
pose, and geometry (ITPG) estimation; 


and (3) refined type, pose, and geometry 
(RTPG) calculation. The image segmen- 
tation module files all the tools in an 
image and isolates them from the back- 
ground. For this, the system uses edge-de- 
tection and thresholding to find the pix- 
els that are part of a tool. After the pixels 
are identified, nearby pixels are grouped 
into blobs. These blobs represent the po- 
tential tools in the image and are the 
product of the segmentation algorithm. 

The second module uses matched fil- 
tering (or template matching) . This ap- 


proach is used for condensing synthetic 
images using an image subspace that 
captu res key i n fo rmati o n . T h ree d egrees 
of orientation, three degreesof position, 
and any number of degrees of freedom 
in geometry change are included. To do 
this, a template-matching framework is 
applied. This framework uses an off-line 
system for calculating template images, 
measurement images, and the measure- 
ments of the template images. These re- 
sults are used online to match seg- 
mented tools against the templates. 
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The final module is the RTPG proces- 
sor. Its role is to find the exact states of 
thetoolsgiven initial conditionsprovided 
by the ITPG module. The requirement 
that the initial conditions exist allows this 
module to make use of a local search 
(whereas the ITPG module had global 
scope). To perform the local search, 3D 
model matching is used, where a syn- 
thetic image of the object is created and 
compared to the sensed data. The avail- 
ability of low-cost PC graphics hardware 
allows rapid creation of synthetic images. 
In this approach, a function of orienta- 
tion, distance, and articulation isdefined 


as a metric on the difference between the 
captured image and a synthetic image 
with an object in the given orientation, 
distance, and articulation. The synthetic 
image is created using a model that is 
looked up in an object-model database. 

A composable software architecture is 
used for implementation. Video is first 
preprocessed to remove sensor anom- 
alies (like dead pixels), and then is 
processed sequentially by a prioritized 
list of tracker-identifiers. 

This work was done by James English, 
Chu-Yin Chang, and Neil Tardella of En- 
ergid Technologies for Johnson Space Center. 


Further information is contained in a TSP 
(see page 1). 

In accordance with Public Law 96-517, 
the contractor has elected to retain title to this 
invention. Inquiries concerning rights for its 
commercial use should be addressed to: 
Energid Technologies 
124 Mount Auburn Street 
Suite200 North 
Cambridge, MA 02138 
Phone No.: (617) 401-7090 
Toll Free No.: (888) 547-4100 
Refer to M SC-23947-1, volume and num- 
ber of this NASA Tech Briefs /ssue, and the 
page number. 


®) Method for Implementing 0 ptical Phase Adjustment 

Goddard SpaceFlight Center, Greenbelt, Maryland 


A method has been developed to me- 
chanically implement the optical phase 
shift by adjusting the polarization of the 
pump and probe beams in an NMOR 
(nonlinear magneto-optical rotation) 
magnetometer as the proper phase shift 
is necessary to induce self-oscillation. 
This innovation consists of mounting 
the pump beam on a ring that sur- 
rounds the atomic vapor sample. The 
propagation of the probe beam is per- 


pendicular to that of the pump beam. 
The probe beam can be considered as 
defining the axis of a cylinder, while the 
pump beam is directed radially. The 
magnetic field to be measured defines a 
third vector, but it is also taken to lie 
along the cylinder axis. Both the pump 
and probe beams are polarized such 
that their electric field vectors are sub- 
stantially perpendicular to the magnet 
field. By rotation of the ring supporting 


the pump beam, itsdirection can be var- 
ied relative to the plane defined by the 
probe electric field and the magnetic 
field to be measured. 

This work was done by David C. H ovdeof 
Southwest Sciences and Eric Corsini of the 
University of California, Berkeley, for God- 
dard SpaceFlight Center. For further informa- 
tion, contact the Goddard Innovative Part- 
nerships Office at (301) 286-5810. 

GSC-15608-1 


@ Visual SLAM Using Variance Grid Maps 

This algorithm is suitable for real-time navigation on irregular terrain. 

N ASA's Jet Propulsion Laboratory, Pasadena, California 


An algorithm denoted "Gamma- 
SLAM " performsfurther processing, in 
real time, of preprocessed digitized im- 
ages acquired by a stereoscopic pair of 
electronic cameras aboard an off-road 
robotic ground vehicle to build accu- 
rate maps of the terrain and determine 
the location of the vehicle with respect 
to the maps. Part of the name of the al- 
gorithm reflects the fact that the 
process of building the maps and de- 
termining the location with respect to 
them is denoted "simultaneous local- 
ization and mapping" (SLAM). Most 
prior real-time SLAM algorithms have 
been limited in applicability to (1) sys- 
tems equipped with scanning laser 
range finders as the primary sensors in 
(2) indoor environments (or relatively 
simply structured outdoor environ- 
ments). The few prior vision-based 


SLAM algorithms have been feature- 
based and not suitable for real-time ap- 
plications and, hence, not suitable for 
autonomous navigation on irregularly 
structured terrain. 

The Gamma-SLAM algorithm incor- 
porates two key innovations: 

•Visual odometry (in contradistinction 
to wheel odometry) is used to estimate 
the motion of the vehicle. 

• An elevation variance map (in contradis- 
tinction to an occupancy or an elevation 
map) is used to represent the terrain. 
The Gamma-SLAM algorithm makes 
use of a Rao-Blackwellized particle filter 
(RBPF) from Bayesian estimation the- 
ory for maintaining a distribution over 
poses and maps. The core idea of the 
RBPF approach is that the SLAM prob- 
lem can be factored into two parts: (1) 
finding the distribution over robot tra- 


jectories, and (2) finding the map con- 
ditioned on any given trajectory. The 
factorization involves the use of a parti- 
cle filter in which each particle encodes 
both a possible trajectory and a map 
conditioned on that trajectory. The base 
estimate of the trajectory is derived 
from visual odometry, and the map con- 
ditioned on that trajectory is a Cartesian 
grid of elevation variances. In compari- 
son with traditional occupancy or eleva- 
tion grid maps, the grid elevation vari- 
ance maps are much better for 
representing the structure of vegetated 
or rocky terrain. 

This work was done by Andrew B. 
Howard of Caltech and Tim K. M arks of the 
University of California San Diego for 
NASA's Jet Propulsion Laboratory. For more 
information, contact iaoffice@jpl.nasa.gov. 
NPO-46114 
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